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FOREWORD 


This  report  documents  work  performed  under  Task  I  of  the  Case-Liner- 
Bond  Study,  authorized  by  Contract  F04611-72-C-0009.  Task  I  involves  the 
calculations  of  typical  case-bond  stress  distributions  in  cylindrical 
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United  States  Air  Force.  The  Air  Force  Project  officer  for  this  work 
is  Mr.  Norman  D.  Walker. 

Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
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exchange  and  stimulation  cf  ideas. 
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ABSTRACT 


Numerical  stress  solutions  were  obtained  for  finite-length  cylinders 
bonded  to  a  flexible  case  and  subjected  to  thermal  shrinkage,  axial 
acceleration,  and  internal  pressure  loading.  Finite-element  models  con¬ 
tained  insulator,  flap,  case-bond  liner,  and  propellant  layers  typical  of 
solid  propellant  rocket  motors.  Grid  structure  in  the  region  of  the  flap 
termination  was  highly  refined  to  provide  accurate  esti^-tes  for  the 
maximum  case-bond  stresses. 

The  stress  solutions  indicated  that  flap  length,  flap  modulus,  and 
liner  modulus  do  not  significantly  affect  the  case-bond  stresses.  The 
most  important  parameter  relative  to  the  maximum  stresses  is  the  thick¬ 
ness  of  material  between  the  flap-insulator  bondline  and  the  liner-pro¬ 
pellant  bondline.  Parametric  studies  were  performed  over  a  range  of 
web-fractions  and  length-to-diameter  ratios.  These  studies  indicate  that 
the  stress  distribution  adjacent  to  the  end  termination  is  insensitive 
to  the  overall  cylinder  configuration,  and  sensitive  only  to  the  local 
end  termination  geometry.  However,  the  maximum  stresses  at  the  end 
termination  increase  with  web-fraction  and  length-to-diameter  ratio  for 
internal  pressure  and  thermal  shrinkage  loading. 

An  equivalence  was  established  between  the  stress  solutions  for 
internal  pressure  loading  and  thermal  shrinkage  loading.  This  equivalence 
provided  highly  accurate  estimates  for  internal  pressure  stresses  using 
the  results  for  restrained  thermal  shrinkage  loading. 
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NOMENCLATURE 

Radius  to  inside  of  grain 

I 

Radius  to  inside  of  case 
C>linder  length 
Cylinder  diameter 
Case-Grain  stiffness  ratio  ( H 


Eb  \ 


Eq  t 
9  c 
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Bond  stress  gradient  factor  (See  Figure  25) 
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SECTION  I 


INTRODUCTION 


Propellant  grain  survivability  depends  upon  the  structural  integrity 
of  'the  grain  and  all  bonded  interfaces.  Probably  more  structural  failures 
of  motors  have  occurred  at  the  case  bond  than  in  the  propellant  proper. 
This  is  true  for  at  least  two  reasons:  (1)  The  case  bond  is  inherently 
weaker  than  the  propellant,  and  (2)  the  stresses  are  highest  at, the  bond 
lines,  particularly  at  the  bond  terminations.  .  , 


Case  bonds  impose  challenges  to  the  structural  analyst  which  are  some¬ 
what  greater  than  those  imposed  by  the  propellant.  These  challenges 
primarily  result  from  .geometrical  considerations.  Layers  comprising  typical 
case-bond  systems,  including  insulator  and  flap  components,  are  very  thin 
compared  with  the  propellant  web.  While  these  thin  layers  do  not  normally 
influence  bond/propellant  stresses  over  most  of  the  motor,  they  are  important 
at  the  critical  bond  termination  locations. 


Only  limited1 information  has  been  published  which  depicts  typical 
case-bond  stress  distributions  in  rocket  motors.  Experimental  studies  have 
been  accomplished  using  photoelasticity  to  evaluate  stress  distributions 
at  bond  terminations.^  However,  these  studies  did  not  consider  case-bond 
detail  and  were  restricted  to  configurations  involving  fillets  and  grooves 
in  the  propellant  grain  adjacent  to  the,  bond  termination.  Numerical  stress 
analyses  published  to  date  have  considered  the  case  bond. in  varying  degrees. 
Reference  2  contains  case-bond  stress  distributions  in  cylinders  with  flat 
ends;  no  attempt  was  made  to  deal  with  the  actual  termination  geometry  or 
bond  constituents.  In  later  studies,-^  however,  some  consideration  was, 
given  to  the  case-bond  liner  and  bond  termination  configuration?.  Fracture 
mechanics  theory  has  been  applied  to  the  case-bond  problem,  with  the  bond 
termination  considered  as  a  singularity^;  information  derived  from  these 
studies  has  been  used  to  establish  the  length  of  flaps..  Finite-element 
stress  analyses  of  motors  often  consider  the  case,  bond;  however,  stress 
values  obtained  at  the  terminations  are  often  of  questionable  value  because 
of  the  lack  of  adequate  grid  refinement. 


Hercules  Incorporated  and  United  Technology  Center  (UTC)  are  currently 
working  on  RPL-funded  programs  to  develop  improved  methods  for  predicting 
initiation  and  propagation  of  failures  in  case-liner-bond  systems.  The 
Hercules  program  involves  initiation  of  case-bond  failure;  whereas,  the  UTC 
program  primarily  involves  propagation  of  case-bond  failures.  The  Hercules 
program  consists  of  four  tasks  involving:  (1)'  Bondline  Parametric  Studies, 
(2)  Test  Sample  Development  for  Simple  States  of  Stress,  (3)  Failute 
Criteria  for  Bond  Terminations,  and  (4)  Subscale  Motor  Testing.  This  report 
documents  studies  performed  under  Task  1. 
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The  objective  of  Task  1,  Bondline  Parametric  Studies,  is  to  establish 
case-bond  stress  distributions  in  basic  motor  configurations,  with  emphasis 
on  bond  terminations.  The  study  involved  simple  cylinders  containing 
flapped  ends.  Finite-element  models  were  developed  with  highly  refined 
grids  in  the  vicinity  of  the  flap  terminations  so  that  variations  in 
case-bond  constituents  could  be  evaluated.  Loading  conditions  considered 
were  thermal  shrinkage,  internal  pressurization,  and  axial  acceleration. 
Information  obtained  under  Task  1  will  be  used  to  develop  meaningful  case- 
bond  failure  criteria  and  associated  test  sample  configurations,  as  outlined 
in  the  remaining  three  tasks.  However,  because  of  its  applicability  to 
propellant  grain  stress  analysis  in  general,  information  gained  under  Task  1 
is  being  published  as  an  interim  report. 

This  report  first  describes  the  modeling  procedure.  Results  are  then 
presented  for  thermal  shrinkage  loading.  An  equivalence  is  developed  for 
thermal  and  pressure  loading  conditions  such  that  thermal  solutions  can  be 
used  for  internal  pressure  loading  conditions  as  well.  Results  are  sub¬ 
sequently  presented  for  internal  pressure  loading  and  compared  with  pre¬ 
dictions  based  on  the  thermal  loading  solutions.  Finally,  results  are 
presented  for  axial  acceleration  loading. 
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SECTION  II 


MATHEMATICAL  MODELING 

The  parameter  study  was  restricted  to  flat-ended  cylinders,  which 
are  symmetric  about  the  axial  mid-plane,  as  shown  in  Figure  1.  The 
cylindrical  propellant  grains  were  bonded  to  the  case,  with  flaps  at  the 
end  terminations.  Emphasis  was  placed  on  modeling  local  detail  in  the 
case-bond  system,  including  the  insulator,  flap,  and  case-bond  liner. 

Stress  solutions  were  obtained  for  the  three  basic  loading  conditions: 

Thermal  shrinkage,  internal  pressurization,  and  axial  acceleration. 

Stress  solutions  were  obtained  using  a  finite-element  computer  program 
based  on  the  quadrilateral  ring  element,  composed  of  four  adjacent  tri¬ 
angular  elements.  The  program  developed  by  Hercules  is  similar  to  the 
well-known  Rohm  &  Haas  program.  ^  The  Hercules  program  is  reformulated  for 
solution  of  problems  with  Poisson's  ratio  near  0.5.  Since  the  reformula¬ 
tion  variable  is  eliminated  through  partitioning  at  the  element  stiffness 
matrix  level,  solutions  cannot  be  obtained  for  r  ~  0.5,  The  program  capacity 
allows  grid  networks  with  up  to  3,000  nodes.  Double-precisioning  is  used 
in  the  program  to  obtain  accurate  stress  solutions  using  the  IBM  System 
360-370  series  computers.  Normal  stress  oscillations  inherent  in  finite- 
element  solutions  using  quadrilateral  elements  were  minimized  through  grid 
refinement  in  areas  of  high  stress  gradient  and  through  averaging  of  stresses 
in  adjacent  elements. 

The  parameter  study  considered  nine  cylinderical  configurations  con¬ 
sisting  of  three  values  of  length-to-diameter  ratio  (L/D)  and  three  values 
of  web  fraction  (W/b) .  The  grid  network  for  L/D  =  1,  and  W/b  =  0.3  is 
shown  in  Figure  2.  The  basic  cylinder  diameter  is  32  in.  to  the  inside 
of  the  case.  Symmetry  conditions  were  imposed  at  the  axial  mid-plane  such 
that  only  one  half  of  the  cylinder  length  was  included  in  the  model. 

Local  detail  in  the  region  of  the  flap  termination  is  indicated  In 
Figure  3.  To  facilitate  modeling  it  was  necessary  to  specify  dimensions  of 
the  cylinders.  The  insulator  is  0.2-in.  thick,  and  the  flap  is  U.l-in. 
thick.  The  flap  length  (length  of  unbondedness  between  insulator  and  flap 
at  the  end  termination)  is  1.0  in.,  and  the  liner  thickness  is  0.06  in. 
Elements  in  the  insulator  are  0.5-in.  thick  (four  elements  to  total  0.2  in.), 
elements  in  the  flap  are  0.025-in.  thick,  and  elements  in  the  liner  graduated 
from  a  minimum  of  0.01-in.  to  a  maximum  of  0. 02-in.  thick.  Elements  in  the 
propellant  adjacent  to  the  liner  are  graduated  in  thickness  away  from  the 
liner,  with  a  minimum  element  thickness  of  0.02  in.  The  element  structure 
in  the  radial  direction  is  carried  through  in  a  uniform  manner  along  the 
total  cylinder  length.  Element  sizes  in  the  axial  direction  are  graduated 
to  provide  minimum  element  lengths  of  0.02-in.  in  the  region  of  the  flap 
termination.  The  grid  contains  2,420  total  nodes. 
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Grid  networks  for  web  fractions  and  lengths  different  from  W/b  =0.8 
and  L/D  =1.0  are  similar  to  that  shown  in  Figure  2.  In  all  grids,  the 
basic  cylinder  diameter  and  insulator  and  flap  thicknesses  were  kept  the 
same.  The  web  fraction  was  changed  by  altering  the  port  diameter,  and  the 
length-to-diameter  ratio  was  changed  by  altering  the  length.  The  grid 
structure  within  4-in.  radially  and  5-in.  axially  of  the  flap  termination 
was  kept  the  same  in  all  grids  for  the  various  web  fractions  and  length-to- 
diameter  ratios.  The  grids  for  L/D  =3.0  were  derived  from  the  grids  for 
L/D  =1.0  using  the  same  number  of  nodes  and  increasing  the  lengths  of  the 
elements  closest  to  the  axial  mid-plane,  where  stress  gradients  are  minimum 
in  the  axial  direction.  The  grids  for  L/D  =5.0  were  otiived  from  the 
grids  for  L/D  =  3.0  by  adding  additional  elements  near  the  axial  mid-plane, 
to  provide  a  total  of  2,640  nodes.  Variations  in  the  web  fraction  were 
obtained  by  decreasing  radial  spacing  between  nodes  inboard  of  the  12-in. 
radial  line  in  the  grids  for  W/b  =  0.8,  without  changing  the  total  number 
of  nodes.  The  thickest  row  of  elements  near  the  center  port,  0.2-in.  thick, 
occurred  in  the  grid  for  W/b  =0.8. 

Material  property  variations  considered  in  the  parameter  study  (Table 
I)  were  kept  to  a  minimum.  The  propellant  tensile  modulus  was  selected 
as  E  =  200  psi,  which  is  approximately  the  longterm  (rubbery)  tensile 
relaxation  modulus  of  many  solid  propellants.  The  propellant  bulk  modulus 
was  selected  as  K  =  333,000  psi,  which  is  the  lower  bound  of  published 
data  for  most  propellants.  The  tensile/bulk  modulus  combination  provides 
a  value  for  Poisson's  ratio  of  v<  =  0.4999  for  propellant. 

The  case-bond  liner  was  assumed  to  have  a  tensile  modulus  equal  to  one 
half  that  of  propellant  for  some  stress  solutions,  but  equal  to  that  of  the 
propellant  in  most  solutions.  The  bulk  modulus  of  the  liner,  flap,  and 
insulator  was  considered  to  be  the  same  as  for  propellant.  (This  choice 
for  the  flap  and  insulator  is  consistent  with  data  for  NBR/SBR  rubber. 6) 

The  tensile  modulus  of  the  insulator  and  flap  was  considered  to  be  equal 
to  that  of  propellant  in  some  solutions  and  equal  to  the  approximate  long¬ 
term  tensile  relaxation  modulus  of  silica-filled  NBR/SBR  rubber  (1000  psi) 
in  most  solutions.  The  shrinkage  coefficient  of  the  insulator,  flap,  and 
liner  was  considered  to  be  the  same  as  for  propellant.  For  thermal  shrinkage 
and  axial  acceleration  loading,  the  motor  case  was  assumed  to  be  rigid. 

For  internal  pressure  loading,  approximate  fiberglass-reinforced  epoxy 
properties  were  used. 
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MATERIAL  PROPERTIES  USED  IN  PARAMETER  STUDY 
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SECTION  ITT 


RESULTS  FOR  THERMAL  SHRINKAGE  LOADING 

Stress  solutions  for  thermal  shrinkage  loading  were  obtained  with  the 
case  (i.e.,  the  outer  diameter  of  the  insulator)  considered  to  be  thermally 
and  mechanically  rigid.  The  propellant,  liner,  flap,  and  insulator  were 
then  assumed  to  shrink  by  the  amount,  b. 

Results  will  be  first  discussed  for  a  typical  set  of  motor  parameters, 
with  case-bond  stresses  shown  as  a  function  of  axial  position,  and  contour 
plots  for  the  region  of  the  bond  termination.  Effects  of  changes  in  the 
flap,  liner,  and  insulator  moduli  on  the  case-bond  stresses  will  then  be 
described.  Results  will  be  subsequently  discussed  for  variations  in  the 
flap  length.  Case-bond  and  ccnterport  stresses  and  strains  will  then  be 
outlined  for  parametric  variations  in  the  grain  web  fraction  and  length- 
to-diameter  ratio. 

Cylinder  parameters  selected  for  the  "typical"  motor  analysis  are: 

L/D  *  1,  W/b  =0.8,  Einsuiator  =  1000  psi,  Efiap  =  1000  psi,  and 
Ejjn^r  =  200  psi  (i.e.,  same  as  propellant).  Stresses  in  the  cylindrical 
coordinate  directions  for  these  properties  are  shown  in  Figure  4  for  the 
row  of  elements  in  the  case-bond  liner  next  to  the  flap.  As  shown  in 
Figure  3,  the  center  of  this  element  row  is  located  0.105  in.  inboard  of 
the  insulator  (0.005  in.  into  the  liner).  The  stresses  were  made 
dimensionless  by  dividing  by  Eb.  These  dimensionless  stresses  are  strictly 
valid  only  for  conditions  wherein  the  insu 1 ator/f lap  modulus  is  five  times 
as  high  as  the  propellant  modulus,  and  for  the  given  values  of  Poisson's 
ra  t  i  o . 

The  stresses  at  the  flap  termination  are  shewn  in  Figure  4;  they  are 
shown  to  be  much  higher  than  at  the  axial  mid-plane  of  the  cylinder 
(Z  -  1 A  in.,  considering  a  32.0-in.  cylinder  diameter).  Radial  stress 
values  averaged  between  adjacent  elements  are  shown  by  the  circles  to 
indicate  grid  spacing  and  stress  gradient  versus  element  size;  point  values 
are  not  indicated  for  the  other  stresses.  The  dimensionless  radial  stress 
near  the  bond  termination  reaches  a  value  of  17.5,  while  it  is  only  3.8 
at  the  axial  mid-plane.  The  dimensionless  hoop  stress  peaks  at  12.3  near 
the  bond  termination,  while  the  dimensionless  axial  stress  peaks  at  8.3. 

The  hoop  and  axial  stresses  are  60  percent  greater  than  the  radial  stress 
at  the  axial  mid-plane.  The  dimensionless  shear  stress  peaks  at  6.3  near 
the  flap  termination  and  goes  through  zero  at  mid-cylinder,  as  specified 
in  the  boundary  conditions  for  symmetry. 

Stress  values  in  the  immediate  vicinity  of  the  flap  termination  are 
shown  in  Figure  5.  Gradients  are  vie1 1  defined  by  the  finite-element 
model.  The  radial  stress  peaks  0.05  in.  aft  of  the  flap  termination 
(i.e.,  a  distance  approximately  equal  to  one-half  the  flap  thickness). 

The  hoop  stress  is  not.  shown,  but  it  peaks  at  the  same  point  as  the  radia1 
stress.  Til-'1  axial  stress  is  still  increasing  along  the  length  after  the 
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radial  and  hoop  stresses  have  peaked.  The  shear  stress  changes  direction 
at  the  flap  termination,  and  peaks  approximately  0.15  in.  aft.  The 
maximum  principal  stress  peaks  slightly  to  the  right  of  the  radiai  °tress 
(as  a  result  of  the  increasing  shear  stress)  at  a  value  of  18.7;  the 
radial  stress,  as  compared  with  the  shear  stress,  is  by  far  the  largest 
contributor  to  the  maximum  principal  stress.  The  maximum  principal  strain, 
fmax/6,  is  also  shown  in  Figure  5;  this  quantity  follows  very  closely 
the  characteristics  of  the  maximum  principal  stress. 

Lines  of  constant  maximum  principal  stress  in  the  vicinity  of  the  flap 
termination  are  shown  in  Figure  6.  The  discontinuity  at  the  bond  termina¬ 
tion  between  the  flap  and  insulator  is  well  defined  by  the  sharp  stress 
gradient.  Since  the  solution  is  linear,  infinite  stresses  should  be  pre¬ 
dicted  at  the  discontinuity;  however,  the  finite-element  model  provides 
stress  values  consistent  with  the  grid  detail  utilized.  The  stress  values 
at  the  flap-liner  interface  are  accurately  estimated  with  the  chosen  grid 
density,  even  though  values  become  more  in  error  as  the  discontinuity  is 
approached.  The  principal  stress  is  not  necessarily  continuous  across 
material  boundaries.  However,  since  the  radial  and  shear  stresses,  which 
must  be  continuous,  are  the  major  contributor  to  the  maximum  principal 
stress,  the  lines  almost  cross  material  boundaries. 

Lines  of  constant  maximum  principal  strain  are  shown  in  Figure  7.  Un¬ 
like  the  stress  values,  the  strain  values  are  quite  discontinuous  across 
the  flap-liner  boundary.  In  the  flap,  there  is  the  expected  concentration 
at  the  flap-insulator  bond  termination  discontinuity;  however,  there  is 
also  an  unexpected  concentration  at  the  flap-lmer  interface.  This  results 
from  a  reduction  in  the  triaxiality  of  the  stresses,  rather  than  by  any 
stress  concentration. 

A.  MATERIAL  PROPERTY  EFFECTS 

The  first  variation  from  the  basic  solution  involved  softening  of 
the  flap  material  from  E  =  1000  psi  to  E  =  200  psi,  making  it  the  same 
as  propellant.  Overall  effects  of  this  change  were  very  small.  Local 
effects  at  the  flap  termination,  as  shown  in  Figure  8,  were  most  signifi¬ 
cant,  but  still  relatively  small.  Radial  and  maximum  principal  stresses 
showed  almost  no  change  from  the  stiffer  flap  solution.  The  shear  stress 
maximum  increased  approximately  6.3  percent  over  the  stiffer  flap  solution. 
The  largest  effects  were  felt  in  the  maximum  principal  strain  and  the 
axial  stress.  With  the  softer  flap  material,  the  axial  stress  decreased 
significantly  at  the  flap  termination,  as  shown  in  Figure  9,  although 
it  built  up  to  the  same  level  as  for  the  stiff  flap  solution  0.5  in. 
further  aft.  This  local  reduction  in  axial  stress  caused  a  more  uniaxial 
stress  condition  to  occur  near  the  location  of  the  maximum  strain,  thus 
increasing  its  value  from  a  maximum  of  10.8  to  13.1.  There  was  an  in¬ 
significant  change  in  the  local  hoop  stress. 
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Lines  of  constant  maximum  principal  stress  and  strain  for  the  soft- 
flap  solution  are  shown  in  Figures  10  and  11,  respectively.  The  maximum 
principal  stress  values  in  the  liner /propellant  are  almost  identical 
for  stiff  and  soft  flap  solutions.  However,  the  stress  concentration  in 
both  the  flap  and  the  insulator  is  significantly  reduced  near  the  bond 
termination  by  softening  the  flap.  Thus,  softening  of  the  flap  has  a 
significant  bearing  on  the  stresses  in  the  rubber  goods  even  though  the 
effect  on  the  bond  stresses  is  minimal.  Lines  of  constant  strain  do  not 
indicate  the  concentration  effect  in  the  flap  adjacent  to  the  1 iner  evident 
in  the  stiff-flap  solution.  The  strain  contour  lines  tend  to  curve  back 
toward  the  concentration  more  than  the  stress  lines,  indicating  the  in¬ 
creased  triaxiality  of  the  stresses  further  aft  of  the  termination. 

The  next  variation  in  material  properties  involved  a  decrease  in  the 
tensile  modulus  of  the  liner,  reducing  it  f. 'm  200  to  100  psi.  For  this 
solution,  the  tensile  modulus  of  the  flap  was  kept  at  200  .psi,  as  opposed 
to  the  standard  value  of  1000  psi.  Case-bond  stresses  and  strains  adjacent 
to  ^  flap  f'  Ruination  for  the  soft  liner  solution  are  shown  in  Figures 
12  through  15;  values  are  plotted  for  the  row  of  elements  in  the  propellant 
nearest  the  liner  (i.e.,  0.01  in.  inboard  of  the  liner-propellant  inter¬ 
face,  and  0.17  in.  inboard  of  the  flap- insulator  interface).  (See  Figure 
3.)  Lines  of  constant  maximum  principal  stress  and  strain  for  the  soft- 
liner  solution  are  shown  in  Figures  16  and  17,  respectively.  In  general, 
the  soft  liner  (in  addition  to  the  soft  flap)  had  little  effect  on  the 
stresses  and  strains  in  the  propellant  or  at  the  liner-propellant  inter¬ 
face.  The  maximum  radial  stress  was  reduced  near  the  liner-propellant 
interface  by  2  percent,  the  shear  stress  by  5  percent,  the  maximum  principal 
stress  by  2  percent,  and  the  maximum  principal  strain  by  7  percent,  as 
compared  with  the  soft  flap  solution. 

The  contour  lines  of  maximum  principal  stress  shown  in  Figure  16  for 
the  soft  liner  solution  indicate  little  change  in  the  stresses  in  the 
liner  itself,  compared  with  the  soft-flap  solution.  However,  the  contours 
for  the  maximum  principal  strain  in  Figure  17  indicate  that  the  liner  is 
shearing  significantly,  several  flap  thicknesses  to  the  right  of  the  flap 
termination.  The  larger  liner  shear  strain,  however,  has  little  effect 
on  the  stresses  and  strain  in  the  propellant  or  at  the  liner-propellant 
interface. 

The  final  variation  in  material  properties  considered  the  insulator 
to  have  the  same  tensile  modulus  as  the  propellant,  liner,  an’  flap 
(E  =  200  psi).  Stresses  and  strains  in  the  row  of  propellant  elements 
immediately  inboard  of  the  liner  are  shown  in  Figures  12  through  15  for 
this  solution.  The  soft-insulator  solution  provided  the  largest  overall 
decrease  in  bond  stresses  and  strains,  i.e.  a  uniform  11-percent  decrease 
in  all  t* e  stresses  and  strains,  as  compared  with  the  soft-flap  solution. 


Contour  lines  of  maximum-  principal  stress  and  strain  for  the  soft 
insulator  solution  are  shown  in  Figures  18  and  19,  respectively.  The 
stress  contours  indicate  a  significant  reduction  in  the  stress  levels  near 
the  insulator-flop  bond  discontinuity,  as  compared  with  the  other  solutions. 

Both  stress  and  strain  contour  lines  are  continuous  across  material  inter¬ 
faces  since  the  material  layers  are  assumed' to  have  the  same  moduli. 

In  summary,  none  of  the  material  property  variations  provided  a 
significant  change  in  the  stresses  at  the  case-bond  interface  adjacent  to 
the  flap  tep» .nation,  except  perhaps  softening  of  the  insulator.  Softening 
of  the  flap  (and  probably  softening  of  the  liner  considering  a  stiffer 
flap)  reduced  the  axial  stress  locally  at  the  flap  termination  such  that 
the  maximum  propellant  strain  increased  significantly,  without  measurably 
affecting  the  radial/shear/maximum  principal  stresses.  Softening  of  the 
liner  increased  the  shear  strain  (and  hence  the  maximum  principal  strain) 
in  the  liner  for  a  distance  to  the  right  of  the  bond  termination,  but 
affected  the  liner/bond  stresses  very  little. 

Perhaps  the  most  significant  feature  of  the  variable-modulus  solu¬ 
tions  is  the  large  stress  gradients,  both  radially  and  axially,  with  dis¬ 
tance  from  the  discontinuity  at  the  insulator-flap  bond  termination.  The 
maximum  principal  stresses  in  each  row  of  elements  are  shown  in  Figure  20 
as  a  function  of  the  radial  position],  normalized  to  unity  at  the  flap- 
liner  interface,  for  the  soft  flap  solution.  The  stress  varies  according 
i  k 

to  the  ratio  ~  ,  where  K  is  a  constant  and  8  is  the  distance  from  the 
\S 

discontinuity.  This  variation  is  consistent  with  the  fracture  mechanics  theory. 
By  doubling  the  flap  thickness,  ‘the  stress  is  reduced  to  72  percent  of 
that-  for  the  standard  0.1-in.  flap  thickness,  and  by  tripling  the  flap 
thickness,  the  stress  factor  goes  to  60  percent.  Also  shown  in  Figure  20 
are  stiff-flap  stresses  which  have  been  normalized  with  respect  to  the 
stress  at  the  flap-liner  interface  for  the  soft  flap.  Only  two  points 
are  plotted  for  the  stiff  flap  solution  since  these  are  sufficient  to 
show  that  stiffening  of  the  flap  had  no  effect  m  the  stresses  inside  the 
propellant.  Corresponding  normalized  stresses  for  the  soft  liner  solution 
indicate  a  similarly  insignificant  effect  on  the  stress  level  inside  the 
propellant.  Thus,  the  propellant  stresses  are  essentially  a  function  of 
the  radial  distance  from  the  flap-insulator  bond  discontinuity  only, 
regardless  of  the  material  layers  (within  the  range  considered)  between 
the  propellant  and  the  discontinuity.  Case-bond  stresses  can  be  reduced 
by  either  thickening  the  flap  or  the  liner.  A  similar  conclusion  cannot 
be  reached  for  the  maximum  principal  strain  since  it  was  shown  to  be 
sensitive  to  the  flap  modulus.  The  stress  behavior,  however,  provides  a 
valuable  procedure  for  estimating  the  effects  of  flap/liner  thickness  on 
the  bond  stresses  using  stress  solutions  fbr  a  single  flap  thickness. 
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R.  FIAP  LENGTH  EFFECTS 

Adjustment  of  the  flap  length  is  a  common  technique  for  minimizing 
propellant/case-bond  stresses.  Information  relating  to  flap  length  effects 
is  available  from  several  sources.  A  simple  procedure  for  sizing  flaps 
was  outlined  by  Lockheed  Propulsion  Company^,  wherein  grain  end  unbonding 
was  observed  in  motors  without  release  flaps.  An  empirical  curve  has  been 
derived  from  motors  with  web  fractions  near  50  percent  and  length-to- 
diameter  ratios  between  1.0  and  3.0.  Based  on  these  data,  a  flap  with 
an  approximate  length  of  1.0  in.  is  desirable  for  a  32- in.  diameter 
cylinder . 

Information  concerning  flap  length  effects  has  been  derived  using 
fracture  mechanics  principles.^  In  fracture  mechanics  solutions,  a 
singularity  is  assumed  at  the  point  of  initial  unbonding.  The  unbond  is 
assumed  to  propagate  when  the  strain  energy  in  the  loaded  motor,  released 
by  a  small  increase  in  the  unbond,  exceeds  the  energy  required  to  create 
a  new  bond  surface  of  the  same  area.  In  general,  as  the  amount  of  un¬ 
bonding  increases  under  thermal  loading,  the  overall  stresses  and  strains 
in  the  motor  decrease,  and  hence,  the  strain  energy  decreases.  However, 
as  shown  in  Reference  4,  the  motor  energy  release  rate  is  a  complex  function 
of  the  motor  web  fraction,  length-to-diameter  ratio,  and  length  of  unbond. 

The  primary  objective  of  the  study  of  flap  length  reported  here  is 
to  evaluate  the  effect  of  varying  flap  length  on  the  local  case-bond 
stress  distribution  at  the  flap  termination.  The  study  is  somewhat 
limited,  with  flap  lengths  of  0.1,  0.5,  and  1.0  in.  considered  for  a 
cylinder  L/D  of  1.0  and  W/b  of  0.8.  Finite-element  models  for  the  two 
shorter  flap  conditions  were  very  similar  to  that  for  the  1.0-in.  flap 
model;  special  attention  was  given  to  keeping  the  grid  sizes  near  the  flap 
termination  identical  in  all  three  models. 

Case-bond  radial  and  shear  stresses  adjacent  to  the  flap  termination 
are  shown  in  Figures  21  and  22,  respectively,  foe  the  three  different 
flap  lengths.  Stresses  are  reported  for  the  row  of  elements  in  the 
liner,  closest  to  the  flap.  The  radial  stresses  for  the  0.1-  and  0.5-in. 
flaps  peak  at  the  same  level;  whereas,  the  peak  for  the  1.0-in.  flap  is 
3  percent  lower.  The  shear  stresses  for  the  0.5-  and  1.0-in.  flaps  peak 
(negative  peak)  at  values  significantly  less  (13  percent)  than  for  the 
0.1-in.  flap.  Overall,  however,  it  is  surprising  that  the  stresses  for 
a  flan  only  as  long  as  its  thickness  are  as  close  as  they  are  to  those  for 
the  longer  flaps. 

The  largest  value  of  the  maximum  principal  stress  in  the  row  of 
elements  in  the  liner  closest  to  the  flap  is  plotted  in  Figure  23  as  a 
function  of  the  normalized  flap  length.  For  the  0.1-in.  flap,  the  largest 
stress  value  is  only  6  percent  greater  than  for  the  l.G-in.  flap.  Thus, 
changes  in  the  bond  stresses  are  relatively  insensitive  to  changes  in 
flap  length  over  the  limited  range  studied  and  flap  thickness  has  a  much 
greater  effect  on  the  case-bond  stresses  at  the  flap  termination  than 
does  flap  length. 
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The  strain  energy  in  the  finite-element  undel  .  shown  in  Firure  24 
as  a  function  of  the  normalized  flap  length.  three  values  calculated 

provide  a  straight  line,  indicating  that  the  ene ■•*>>•  release  rare  is  con¬ 
stant.  This  is  consistent  with  local  stresses  be  ’  >  al-~r.t  independent 
of  flap  length;  differences  between  energy  and  local  stress  criteria 
can  probably  be  explained  by  the  limited  number  of  solutions  used  to 
define  the  energy-versus-flap-length  curve.  Considering  g  -  1.0  in., 
flap  length  effects  on  the  local  case-bond  stresses  for  cylinders  with 
smaller  web  fractions  and  larger  L/P  ratios  are  expected  to  be  tpss  than 
those  for  the  stubby  cylinder  which  was  studied,  based  on  energ'r  release 
rate  curves  contained  in  References  2  and  4. 

C.  WEB  FRACTION  AND  LENGTH /DIAMETER  (L/D)  STUDIES 

Overall  case-bond  stress  distributions  in  cylinders  have  been  shewn  to 
be  sensitive  to  web  fraction  and  length-to-diameter  ratio.  (See  Reference  2 
for  exutrple.)  It  is  important,  therefore,  that  the  current  study  consider 
as  an  objective  the  evaluation  of  case  bond  stress  distributions  for  varia¬ 
tions  in  W/b  and  L/D.  A  secondary  objective  is  to  provide  design  charts 
which  car  be  used  to  estimate  maximum  values  of  case-bond  stresses  at  flap 
terminations . 

The  study  considered  nine  cylinder  configurations:  L/D  =1,3,  and 
5,  for  W/b  =  0.5,  0.67,  and  0.8.  The  cylinder  diameter  was  held  constant 
at  32  in.  for  these  studies;  thus,  L/D  was  varied  by  changing  L  and  W/b 
was  varied  by  changing  propellant  bore  diameter.  The  thickness  and  modulus 
of  the  flap  and  ir.iuiator,  the  length  of  the  flap,  and  the  modulus  of  the 
propellant  were  kept  the  same  in  all  nine  solutions.  Modulus  values  used 
in  the  study  are  indicated  in  Table  1.  The  insulator  thickness  of  0.2  in. 
is  the  same  as  for  all  previous  solutions,  and  the  flap  thickness  of  0.1 
in.  and  the  length  of  1.0  in.  is  also  the  same  as  that  used  in  previous 
solutions . 

Data  shown  in  Figure  20  indicate  that  the  propellant  stress  adjacent 
to  fhe  flap  termination  is  not  strongly  affected  by  the  modulus  of  the 
flap  or  liner.  Thus,  the  stresses  obtained  in  a  computer  solution  for  a 
given  flap/liner  thickness  can  he  made  to  apply  to  conditions  in  which 
the  flap  and  liner  thickness  vary.  Bond  normal  stresses  are  presented 

parametrically  as  ~,  where  Ho-  is  a  function  of  the  radial  distance,  t, 

of  the  prope llant-to-_liner  bondine  from  the  flap-to-insulator  bondline;  for 
a  liner-bond  system,  t  =  tflap  +  clir  r  wben  the  critical  failure  location 
is  between  the  propellant  and  liner.  The  value  for  is  presented  in 
Figure  25  as  a  function  of  the  normalized  thickness,  I/b.  The  value 
was  normalized  such  that  if  would  be  unity  at  a  distance  of  0.x,,..  in. 
from  the  flap-insulator  bondline,  where  stresses  are  obtained  conveniently 
from  che  finite-element  solutions.  The  value  for  HT  may  be  a  function  of 
W/b  and  L/D;  the  value  shown  in  Figure  25  is  for  L/D  =1.0  and  W/h  =  0.8. 
Figure  26  contains  a  comparison  of  Ho-  for  the  extreme  values  of  cylinder 
parameters  studied  (i.e.,  L/D  =  3,  W/h  =  0.8,  and  L/D  =  5,  W/b  =  0.5)  . 


There  is  a  negligible  difference  in  H0  for  these  extreme  conditions,  indi¬ 
cating  that  the  stress  distribution  adjacent  to  the  flap  termination  is  a 
function  of  the  local  flap-insulator  geometi'y  and  not  the  overall  cylinder 
geometry.  Furth*.  rmore,  the  radial  stress  gradient  near  the  flap  termina¬ 
tion,  as  shown  in  Figure  26,  is  so  close  to  that  for  the  maximum  principal 
stress  that  in  Figure  25  applies  to  both  stress  components,  although 
it  was  obtained  from  the  maximum  principal  stress.  As  shown  in  Figure  26, 
Ho-  for  normal  stresses  does  not  apply  to  the  rrz  shear  stress.  The  peak 
value  of  the  shear  stress  does  not  decrease  as  rapidly  with  distance  from 
the  discontinuity  as  do  the  normal  stresses. 

A  relationship  for  the  strain  components  as  a  function  of  flap  thick¬ 
ness  cannot  be  obtained  with  confidence  from  a  solution  for  one  flap 
thickness  since  the  flap  modulus  significant!-  affects  the  strain  magni¬ 
tude,  as  shown  in  Figures  5  and  8.  Results  are  provided,  however,  for  the 
basic  flap  thickness  to  indicate  variations  with  cylinder  web  fraction  and 
length . 


The  dimensionless  radial  stress  adjacent  to  the  flap  termination  is 
shewn  in  Tigure  27  as  a  function  of  L/D  for  the  three  web  fractions. 

Values  change  significantly  between  L/D  =1.0  and  3.0,  and  little  there¬ 
after;  the  greatest  effect  occurs  for  the  largest  wee  fraction.  Also 
indicated  in  Figure  27  is  the  ratio  of  the  shear  stress  rrz  to  the  radial 
stress  at  the  location  of  the  maximum  radial  stress  for  the  standard  flap/ 
liner  thickness,  "t  =  0.105  in.  The  stress  ratio  changes  very  little  with 
web  fraction  and  length,  indicating  that  the  local  stress  distribution 
(sh<.„n  in  detail  earlier)  for  L/D  =1.0  and  W/b  =  0.8  is  characteristic 
of  other  cylinder  configurations  as  well. 

To  illustrate  the  use  of  parametric  information  in  Figure  27  for  design 
purposes,  consider  a  motor  with  the  following  dimensions: 


Item 

Case  Diameter,  D 
Flap  Thickness,  r. 
Liner  Thickness,  t, 
L/D 
W/b 


Value 

12  in. 

0.08  in. 
O.05  in. 

2 

0.75 


t.  .  /b 

insulator 

E  /E 

insulator  propellant 


propellant 

v'insulator 

vflap 


v  liner 


0.2/16  =  0.0125 
5 

0.4999 

0.4995 

0.4999 

0.4999 
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The  interpolated  value  ot  I— — g^jbond  ^rora  Flgure  27  is  approximately  24.0. 
The  normalized  sum  of  the  liner  and  flap  thickness  is, 

t/b  =  (0.08  +  0 .05) /I2  =  0.0108  ; 

thus,  H(7  from  Figure  25  is  0.79.  (The  critical  case-bond  stress  location 
is  assumed  tc  be  between  the  liner  and  propellant,  not  between  the  flap 
and  liner.)  Assuming, 

E  =  500  psi 

6=  a  AT  -  0.004  in. /in., 


then 


cr  =  (0.19)  (500)  (0.004)  (24)  =  37.9  psi 

The  case-bond  snear  stress  at  this  location  is  approximately 

t  .  -  (0.19)  (37.9)  =  7.2  psi 
r  o 

The  dimensionless  maximum  principal  stress  adjacent  to  the  flap 
termination  is  shown  in  Figure  28  and  the  maximum  principal  strain  is 
shown  in  Figure  29,  The  strain  is  located  at  T/b  =  0.0066  only,  and 
varies  with  v:eb  traction  and  length  similar  to  the  stresses.  Also  shewn 
in  Figure  29  are  the  ratios  of  the  other  two  principal  stresses  to  the 
maximum  principal  stress  at  the  location  of  the  maximum  principal  strain 
(same  location  as  the  maximum  principal  stress  for  T/b  =  0.0066).  The 
hoop  stress  is  approximately  65  percent  of  the  maximum  principal  stress, 
and  the  smallest  principal  stress  is  approximately  20  percent  of  the 
maximum  principal  stress.  There  is  only  a  small  variation  in  the  stress 
ratios  with  L/D  and  W/b, 

For  completeness,  the  centerport  hoop  stress  and  strain  at  the  mid¬ 
cylinder  is  presented  in  Figures  30  and  31,  respectively.  These  curves 
are  almost  identical  to  those  presented  ir.  Reference  8,  wherein  finite 
difference  methods  were  used  foi  the  calculations.  The  strain  values 
are  for  the  stress-producing  strain,  «“«•*■  6,  as  opposed  to  the  total 
strain,  c=  ur/r.  The  centerport  stress  and  strain  do  not  level  off  as 
quickly  with  increasing  L/D  as  do  the  bond  termination  stresses  and 
strains.  The  ratio  of  the  &,.Ial  stress  (a  principal  stress)  to  the  hoop 
stress  at  the  mid-cylinder  is  indicated  in  Figure  31.  This  stress  ratio 
changes  dramatically  with  L/D  and  W/b,  and  is  approximately  0.5  (i.e., 
2:1  biaxial)  tor  L/D,  greater  than  3.0,  ana  approaches  zero  (uniaxial) 
for  L/D  equal  to  1.0, 
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The  ratio  of  the  maximum  principal  stress  at  the  flap  termination  to 
the  maximum  centerport  hoop  stress  is  shown  in  Figure  32.  As  the  web 
fraction  and  length-to-diameter  ratio  decrease,  the  bond  termination 
stress  becomes  much  larger  chan  the  centerport  hoop  stress.  For  high  L/D 
cylinders  with  large  web  fractions,  the  centerport  hoop  stress  becomes 
larger  than  the  bond  termination  stress,  although  the  bond  termination 
stress  can  be  adjusted  appreciably  by  the  flap/liner  thickness  (Ho-)  . 

The  case-bond  radial  stress  distribution  as  a.  function  of  motor  length 
varies  greatly  with  L/D,  as  indicated  in  Figure  33.  For  L/D  =1.0  and 
W/b  =  0.08,  the  radial  stress  at  the  mid-cylinder  is  muchness  than  at  the 
flap  termination.  However,  for  L/D  =  3.0,  the  radial  stress  increases 
significantly  near  che  mid-cylinder,  and  a  distinct  saddle  region  of 
minimum  stress  occurs  to  the  right  of  the  flap  termination.  The  ratio 
of  the  maximum  case-bond  radial  stress  at  the  flap  termination  to  the  case- 
bond  radial  stress  at  the  mid-cylinder  is  shown  in  Figure  34  as  a  function 
of  L/D  and  W/b.  The  trend  for  this  ratio  is  similar,  to  that  for  the  flap- 
centerport  hoop  stress  ratio;  however,  the  radial  bond  stress  is  less  than  1 
the  centerport  hoop  stress  at  the  mid -cylinder  location.  The  ratio  of  the 
case-bond  radial  stress  at  the  saddle  to  that  at  the  mid-cylinder  is  shown 
in  Figure  35.  The  saddle  deepens  with  increasing  L/D  and  W/b. 

It  was  shown  in  Figure  4  that  the  hoop  and  axial  stresses  along  the 
case-bond  line  are  greater  than  the  radial  stress,  except  near  the  flap 
termination.  This  may  be  somewhat  surprising  at  first  glance.  However, 
the  plane  strain  cylinder  solution  for  thermal  loading  predicts  that 


for  v  =  C.5  and  a  rigid  case.  Thus,  as  the  web  fraction  decreases,  tne 
hoop  stress  at  the  bond  becomes  larger  than  the  radial  stress.  The 
finite-element  results  for  the  ratio  of  the  hoop  stress  to  the  radial 
stress  at  the  mid-cylinder  bond  location  are  given  in  Figure  36.  The 
dependence  upon  web  fraction  for  large  values  of  L/D  (approaching  plane 
strain)  follows  the  plane  strain  predictions.  In  addition,  the  hoop 
stress  increases  in  magnitude  over  the  radial  stress  'as  L/D  decreases. 
The  axial  stress  at  the  mid-cylinder  is  equal  (within  1-percent  error) 
to  the  hoop  stress. 

The  average  case-bond  radial  stress  is  shown  in  Figure  37  for  the 
various  cylinder  parameters.  The  distribution  indicates  an  increasing, 
average  load  cn  the  case  with  increasing  L/D  and  W/b. 
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SECTION  IV 


THERMAL- PRESSURE  LOADING  EQUIVALENCE 

Much  similarity  exists  between  the  stress  solutions  for  continuums 
subjected  to  restrained  thermal  shrinkage  and  pressure  loadings.  (See 
Reference  9,  for  example.)  This  similarity  can  be  used  to  advantage  in 
both  stress  and  strength  analysis  of  solid  propellant  rocket  motors.  In 
particular,  the  number  of  stress  solutions  required  to  depict  typical 
case-bond  stress  distributions  in  tWo-layered,  finite-length  cylinders 
can  be. reduced  in  the  current  parameter  study  by  taking  advantage  ol  this 
equivalence. 

A  pictorial  presentation  of  propellant  (and  case  bond)  stresses  and 
strains  in  cylindrical  rocket  motors  with  free  ends,  subjected  to  thermal 
shrinkage  and  internal  pressure  loading,  is  given  .in  Figure  38.  During 
thermal  shrinkage  loading,  propellant  normal  stresses  f<r),  shear  stresses 
(r),  and  stress-producing  strains  (?)  occur  because  the  propellant  has  a 
higher  coefficient  of  thermal  expansion  than  the  case. 

The  stresses,  strains,  and  deformations  for  thermal  leading  can  be 
separated  into  those  resulting  from  two  situations,  free  shrinkage  and 
case  expansion,  as  indicated  in  Figure  38-a.  Assuming  that  the  propel¬ 
lant  is  disconnected  from  the  case  and  allowed  to  uniformly  shrink, 
zero  normal  and  shear  stresses  will  result.  By  definition,  the  stress- 
producing  strains  (?)  will  also  be  zero;  the  total  strains  (e)  will  then 
be  equal  to  the  negative  of  the  linear  shrinkage  coefficient  (<5p)  .  The 
case  will  also  shrink,  but  normally  much1 less  than  the  propellant.  When 
the  case  and  propellant  are  bonded  togethu. ,  interface  stresses  will 
occur  at  the  case-bond  line.  The  case  deformation  resulting  from  these 
stresses  Is  usually  small,  but  is1 dependent  upon  the  relative  grain-to- 
case  stiffness  ratio  in.  the  coupled  system,  as  are  the  interface  stresses. 
For  convenience  of  presentation,  consider  that  the  case  deformation  due 
to  the  interface  stresses  is  known  and  included  in  the  free  shrinkage 
solution.  The  propellant  stresses  for  restrained  shrinkage  loading  can 
then  be  viewed  as  those  which  occur  when  the  outer  boundafy  of  the  shrunk 
propellant  cylinder  is  pulled  out  to  make  contact  with  the  deformed  case 
(i.e.,  case  expansion  loading),  as  indicated  in  the  right  side  of 
Figure  38-a.  Since  the  propellant  grain  under  free  shrinkage  loading  ^as 
assumed  to  uniformly  shrink  such  that  all  total  normal  strains  are 
equal  to  "-5p",  and  the  case  deformation  is  usually  small  due  to 
the  interface  stresses,  the  hoop  strain  (f q)  will  be  approximately  equal  to 
the  axial  strain  ((%)  all  along  the  case-bond  line  for  case  expansion 
loading.  Strength  analysis  is  most  logically  carried  out  in  terms  of 
stress  producing  strains,  f,  as  opposed  to  total  strains,?. 
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During  uniform  internal  pressure  loading,  propellant  stresses  (other 
than  hydrostatic  pressure  stresses)  and  strains  occur  because  the  case 
expands  and  the  grain  follows  to  take  up  the  increased  volume  inside  the 
deformed  chamber.  The  stresses,  strains,  and  deformations  for  uniform 
internal  pressure  loading  can  be  separated  into  those  resulting  from  two 
situations,  hydrostatic  pressure  and  case  expansion,  as  indicated  in 
Figure  38-b.  Assuming  that  the  propellant  is  disconnected  from  the  case 
and  subjected  to  uniform  hydrostatic  pressure,  all  normal  stresses  will 
be  equal  to  -P  and  all  shear  stresses  will  be  zero;  the  normal  strains 
P 


will  be  equal  to 


3K‘ 


Consider  now  a  strain  definition  corresponding  to 


that  for  the  stress-producing  strain,  7,  under  thermal  loading.  In  the 
pressure  loading  situation,  this  "stress-producing  strain,"  7,  corresponds 
to  that  required  to  produce  normal  stresses  different  from  -P;  thus,  7 
by  definition  is  zero  for  hydrostatic  pressure  loading. 


For  convenience  of  presentation,  consider  that  the  case  deformation 
due  to  interface  stresses  and  end  forces  is  known  for  internal  pressure 
loading,  even  though  it  is  a  function  of  the  coupled  grain/case  system. 

The  propellant  stresses  considering  the  propellant  and  case  are  bonded 
together  can  then  be  viewed  as  those  for  hydrostatic  pressure  loading, 
plus  those  which  occur  when  the  outer  boundary  of  the  compressed  propellant 
cylinder  is  pulled  out  to  make  contact  with  the  deformed  case.  The  case 
strains  for  an  empty  metal  chamber  will  be  such  that  (■  7  „  ^  f,(>  consider¬ 
ing  a  filament-wound  glass,  *£  will  be  typically  as  large  as  f q.  Con¬ 
sidering  a  case  designed  to  produce  a  hoop  strain  of  0.01  in. /in.  at 
1000  psi  internal  pressure,  and  a  propellant  with  a  bulk  modulus  of  333,000 
psi,  the  propellant  strain  due  to  uniform  hydrostatic  compression  will  be 

£  =  =  -0.003.  in. /in.,  which  is  10  percent  of  the  empty  case  hoop 

strain.  The  empty  case  expansion  component  of  loading  can  usually  be 
expected  to  dominate  the  propellant  compressibility  component  of  loading. 
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SECTION  V 


LOAD  SHARING  OF  CASE  AND  GRAIN 

The  interface  stresses  between  the  case  and  grain  are  a  function  of 
the  relative  grain-to-case  stiffness.  A  procedure  for  estimating  the  case 
deformation  which  results  from  the  interface  stresses  is  derived  in  the 
following  paragraphs.  Pressure  loading,  which  produces  the  most  complex 
case-grain  coupling,  is  discussed  in  detail;  then,  case-grain  coupling 
is  summarized  for  thermal  loading. 

A.  INTERNAL  PRESSURE  LOADING 

Consider  a  two-layered  cylinder  subjected  to  internal  pressure  loading 
The  case  hoop  and  axial  stresses  can  be  approximately  expressed  as: 


a 


e 


(p  -  <r, 


bond 


<r 


r 


bond 


The  term  ,,crbond"  rePresents  the  average  internal  pressure  load  carried 
by  the  grain,  andTboncj  relates  to  the  shear  stress  at  the  interface. 

The  pressure  load  carried  by  the  grain  is  highest  near  the  end  terminations 
where  there  is  a  sharp  stress  concentration,  and  near  the  axial  mid-plane 
where  gross  end  relief  effects  diminish.  The  shear  stress  is  highest 
near  the  end  terminations  and  changes  sign  at  the  axial  mid-plane  such  that 
the  integral  of  the  shear  is  zero  over  the  total  length.  Shear  is 
secondary  compared  with  the  radial  stress  effect,  and  will  be  ignored  in 
further  considerations. 


For  a  membrane  case,  the  orthotropic  stress-strain  relationships  are: 


”  cse  °e  +  cez  °z 


'z  *  CZ9  °i5  +  czz 


Substitution  of  the  above  stress-load  relationships  (minus  the  shear) 
into  *-he  stress-strain  relations  provides: 


f9  =  C09  '  °bond^  tc  +  C9Z  (ztc  J 

?Z  =  CZ8  (P  '  flrbond)  t  +  GZZ  (fr) 
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Upon  regrouping  these  expressions,  the  following  are  obtained: 


e  =  e  -  r  —  it 
9  0  00  t  bond 

c 


£  =  (  -  c  —  <T 

Z  6  60  t  bond 

c 


where  eQ  and  t  are  the  case  strains  in  the  absence  of  the  grain. 

U  L, 

The  average  interface  radial  stress  carried  by  the  grain,  °i,ond> 
is  a  function  of  both  the  radial  and  axial  case  strains.  Consider  functions 


®  a"d  (%) 


which  are  the  average  dimensionless  bond  radial  stresses  for  arbitrary 
(case  expansion  loading)  hoop  and  axial  case  strains.  Further,  recognize 
1  V  QZ 

CrtQ  =  7“  anc^  C  a  =  -  — —  .  Then  the  hoop  and  axial  case  strains  are 
Bo  Eg  ZB  Eg 

related  as: 


=  '  e9  tc  LUe 


f(it)  ^ +  (^)  <z. 


£  =  f  +  V  — — 

Z  Z  0Z  E0  tc 


-(i'e)  v(.r0  <z] 


These  equations  can  be  solved  for  £g  and  as: 


»  f9-  (,V!Z _ 

-  V.  H  »■«)  (1  +  11  °rB*  +  V  “2  CTrZ  CTrO 


TZ  +  (v0Z  H  ar0)  *  0 


1  -  "ez  H  CrZ 
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where: 


H  = 


Bb  (  °r  \ 

=9  'c  ’  r0  "  lEfe 

case ' 

avg 


and  <t 
rZ 


.  U 


CT 


case^ 


avg 


B.  THERMAL  SHRINKAGE  LOADING 

Case  axial  and  hoop  stresses  for  thermal  shrinkage  loading  can  be 
approximately  expressed  as: 


<-  *W 


b_ 

t 

c 


o 


The  stress  producing  case  strains  arc  related  to  the  average  bond  radial 
stress  as: 


e 

0 


°bond^ 


b_ 

t 

c 


€ 

Z 


0Z 


^bond^  t 

c 


U9Z  f0 


(<rv\ 

The  average  bond  radial  stress ,1  rtl  ,  is  presented  in  Figure  37  for 

'  Eo/avg 

variations  in  L/D,  W/b,  E,  and  5  of  a  cylinder  with  a  rigid  case.  Thus, 


and: 


UQZ  f0 
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The  stresses  for  thermal  shrinkage  loading,  considering  a  flexible 
case,  are  then  approximately  equal  to  those  for  a  rigid  case  plus  those 
for  case  expansion  loading  of  the  amount  estimated  above.  Since  the  case 
strains  are  negative,  however,  the  flexible  case  solution  will  provide 
lower  stresses  than  the  rigid  case  solution. 
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SECTION  VI 


RESULTS  FOR  INTERNAL  PRESSURE  LOADING 

A.  DIRECT  INTERNAL  PRESSURE  LOADING  RESULTS 

The  case-bond  stress  distribution  for  a  typical  set  of  cylinder 
parameters  was  obtained  for  internal  pressure  loading  applied  to  the  finite- 
element  model.  The  propellant,  liner,  and  insulator  properties  used  are 
those  underlined  in  Table  I.  However,  the  flap  modulus  was  200  psi  and 
the  motor  case  was  not  assumed  to  be  rigid;  a  case  with  thickness,  tc  =  0.1 
in.  (tc/b  =  0.00625),  a  tensile  modulus,  Ec  =  4  x  10“  psi,  and  Poisson's 
ratio,  vc  =  0.1  were  used.  These  case  properties  provide  a  mid-cylinder 
case  hoop  strain  of  1.85  percent  at  500  psi  internal  pressure,  which  is 
typical  of  fiberglass -reinforced  epoxy  cases. 

Case-bond  stresses  in  the  row  of  elements  in  the  case-bond  liner  next 
to  the  flap  in  the  immediate  vicinity  of  the  flap  termination  are  shown 
in  Figure  39  for  P  =  1.0  psi.  This  stress  distribution  is  virtually 
identical  with  the  companion  solution  for  thermal  shrinkage  loading  shown 
in  Figure  5.  Lines  of  constant  maximum  principal  strain  are  shown  in 
Figure  40  for  the  immediate  vicinity  of  the  flap  termination;  gradients 
are  almost  identical  with  those  for  the  corresponding  shrinkage  solution 
in  Figure  11.  Thus,  there  is  little  value  in  presenting  more  of  the 
stress  distributions  for  internal  pressure  loading  since  they  are  similar 
to  those  already  presented  for  thermal  shrinkage  loading. 

B.  CASE  HOOP  EXPANSION  COMPONENT  OF  LOADING 

Stress/strain  values  already  presented  for  thermal  shrinkage  loading 
can  be  interpreted  as  case  expansion  solutions  for  =  fi  -  8  at  the 
case-insulator  interface.  By  obtaining  additional  stress  solutions  for 
a  unit  value  of  fg  of  the  case  and  a  zero  the  stresses  for  arbitrary 
combinations  of  f q  and  can  be  derived.  Key  stresses  for  f q  case 
expansion  loading  of  the  cylinder  with  an  80-percent  web  fraction  are 
presented  in  Table  II.  They  are  compared  with  the  corresponding  stresses 
for  thermal  shrinkage  loads;  however,  the  shrinkage  stresses  were  reduced 
through  division  by  a  factor  of  1.5.  The  1.5  factor  is  based  on  the 
assumption  that  the  cylinder  stresses  for  a  unit  axial  strain  are  equal 
to  Poisson's  ratio  (~0.5)  times  those  for  a  unit  hoop  strain.  For  thermal 
shrinkage  loading,  the  axial  and  hoop  strains  were  assumed  to  be  equal 
at  the  case,  considering  shrinkage  to  be  case  expansion  loading.  Thus, 
stresses  for  hoop  case  expansion  loading  only  should  be  two-thirds  those 
for  shrinkage  loading.  The  comparison  between  the  (reduced)  shrinkage 
values  and  the  hoop  case  expansion  values  is  surprisingly  good,  with 
deviations  generally  less  than  10  percent.  The  largest  deviations  occur 
for  the  centerport  hoop  stress  and  strain  for  L/D  =  1  (14  percent),  and 
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TABLE  II 


1 


for  ratios  of  the  flap  and  saddle  stresses  to  the  mid-cylinder  radial 
stress  at  the  bond  for  L/D  =  5.0  (16 ' percent) .  The  critical  propellant 
and  bond  stresses  for  case  hoop  expansion  loading  can  be  predicted  with 
good  accuracy  from  the  thermal  shrinkage  solutions.  Therefore,  case  hoop 
expansion  solutions  were  not  obtained  for  the  other  two  web  fractions. 

C.  '  PROCEDURE  FOR  PREDICTING  PRESSURE  SOLUTIONS  FROM  THERMAL  SOLUTIONS 


Stresses  and  strains  for  internal  pressure  loading  can  be  predicted 
completely  from  corresponding  values  for. thermal  shrinkage  loading.  The 
procedure  for  predicting  stress  solutions  for  internal  pressure  loading  is 
as  follows: 

(1)  Calculate  axial  and  hoop  strains  for  the  case,  without  the 
grain:  '  t 


case 


case 


E 


0 


cr 

2 


I 


where  <T ^ 


■ —  and  cr 
t  Z 

c 


Pb 

2t 


(2)  Determine  average  case  strains  for  coupled  case  and  grain: 


rx .  „  H  ( -r.\  i  ?  .  h?  i 

ez  3  \i :dj  e  |3  v F§/ 
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3  0z'H\eSJ  J 

-  I  V  H  fe) 
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fZ  + 


ravg  case 


where  H  =  ;  - — —  ,  and  J  is  obtained  from  Figure  37, 
9c  '  'avg 


J 


(3)  Calculate  critical  stresses  and  strains  for  pressure  loading: 


a  +  P 

pressure 


e 

9 


case 


( 

Z 

case 
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P_ 

3K 


hrinkage 


€ 

pressure 


-  €  +  -  f  4. 

Gc»se  3  '  3K 


hrinkage 


Note:  Case  strains  for  the  empty  case  must  apply 
for  pressure,  P. 

This  procedure  will  now  be  demonstrated  for  the  cylinder  with  W/b  =  0.8 
and  L/D  =  3.0,  considering  an  isotropic  case. 

(1)  Empty  case  strains: 


fr  _  Ik  _  (1.0)  (16.05)  _ 

0  t  0.1 

c 


=  160.5 


Pb 

2t 


=  80.25 


-  cr  .  _£ 
9  °  E0 


a 
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4  x  10  4  x  10 


g  (80.25)  =  38.1  x  10 
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^  <r  +!—  <7  = 

e0  e  ez  Z 
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-6 
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(2)  Case  strains  considering  grain  reinfor.  ement : 


H  = 


Eb 


(200)  (16) 


E9  Cc  (4  x  106)  (0.1) 


=  0.0080 


_r 

ES, 


=  11.2 
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1 

£  — 

1  -  (0.1) (0.0298) J 

(38.1  x  10“6)  -(0.0298)06.04  x  10_6) 

0 

case 

1  -  (0.1) (0.0298) j 

(1  +  0.0596)  -f  (2)  (0.1)  (0.0298)2 

=  35.6  x  10‘6 

16.04  x  10~6  +  (2)  (0.0298)  (35.6  x  10~6)  =  -6 

1  -  0.1  (0.0298) 

(3)  Calculate  critical  stresses  and  strains: 


(CTrt  +  P).  =  200  [ (0.667) (35.6  x  10~6)  +  (0.333)(18.2  x  10  6) 

0  bore  l 

+  1  x  10'6  j  (32.6)  =  0. 


(cr  +  P)„  =  200  (30.8  x  10"6)  (33.0)  =  0.203 

max  r  lap 

9)bore  =  (30'8  X  10"6)  (24‘9)  =  °-00077 


A  finite-element  solution,  obtained  directly  for  internal  pressure 
loading,  correlated  with  the  above  prediction  as  follows: 


Stress 

Pressure  Loading 

Parameter 

Result 

(<r9  +  P)bore 

0.191 

^  +  P>£lap 

0.200 

(<9'bore 

0.0074 

Predicted  Value 
from  Thermal  Solution 

Error 

(percent) 

0.201 

+5.2 

0.203 

+1.5 

0.00077 

+3.5 
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The  maximum  error  of  5.2  percent  occurs  in  the  centeiport  hoop  streSo.  In 
all  cases,  the  thermal  prediction  is  high.  The  predicted  case  hoop  strain 
of  35.6  x  10-6  compares  with  pressure  loading  results  of  35.7  x  10"^  at  the 
flap  termination,  36.6  x  10"®  at  the  saddle  (maximum  case  hoop  strain),  and 
34.5  x  10"6  at  the  mid-cylinder;  the  predicted  case  axial  strain  of  13.2 
x  10“^  compares  with  pressure  loading  results  of  16.6  x  ] at  the  flap 
termination,  1C - 1  at  the  saddle,  and  15.0  x  10“^  at  the  mid-cylinder. 

Thus,  the  average  hoop  strain  in  the  case  compares  quite  well  with  the 
prediction,  while  the  predicted  axial  strain  is  too  h_gh .  Omission  of 
the  bond  shear  stress  in  the  grain-case  coupling  equations  is  probably 
to  blame  for  the  high  axial  strain  prediction.  Overall,  the  directly 


calculated  pressure  solution  compares  very  favorably  with  the  prediction 
based  on  thermal  shrinkage  loading. 


SECTION  VII 


RESULTS  FOR  AXIAL  ACCELERATION  LOADING 


Sti'ess  solutions  for  axial  acceleration  loading  were  obtained  with 
the  case  considered  to  be  ligid.  The  propellant,  liner,  flap,  and 
insulator  were  given  the  same  weight  density,  =  0.064  Ib/in.^  and 
subjected  to  a  1.0  g  axial  acceleration.  Stress  solutions  were  obtained 


only  for  £.  ,  _  =  EC1 

insulator  flap 

The  flap  length  was  1.0  in. 


1000  psi  and  E 


liner  ^propellant 


=  200  psi. 


Calculated  radial  and  shear  stresses  for  the  cylinder  with  L/D  =  3.C 
and  W/b  =  0.8  are  shown  in  Figure  41  for  the  row  of  elements  in  the  case- 
bond  liner  next  to  the  flap.  The  center  of  this  element  row  is  located 
0,105  in.  inboard  of  the  insulator  (0.005  in.  into  the  liner).  The  shear 
stress  reaches  an  average  value  of  -0.50  psi  at  the  mid-axial  location 
(Z  =  4£  in.).  Near  the  flap  termination  the  shear  stress  has  a  negative 
peak  of  -0.68  psi  and  '  positive  peak  of  0,23  psi.  The  peak  radial 
stress  at  the  flap  termination,  3.4  psi,  is  much  higher  than  the  shear 
stress.  Both  the  low  value  of  the  peak  shear  stress  and  the  high  value 
of  the  peak  radial  stress,  are  somewhat  surprising.  Currently-used 
handbook  methods  for  approximating  stresses  at  the  bond  termination  assume 
that  the  peak  shear  stress  is  3.0  times  the  average  shear  stress.1®  In 
the  sample  given  here,  the  calculated  concentration  ictor  for  the  peak 
shear  stress  is  only  1.4.  However,  the  peak  radial  stress  is  4.3  times 
the  average  shear  stress. 


Stress  values  in  the  immediate  vicinity  of  the  flap  termination  are 
shown  in  Figure  42.  Ihe  shear  stress  and  maximum  principal  strain  approach 
their  average  values  to  the  right  of  the  flap  termination  rather  than 
approaching  zero,  as  in  the  thermal  solution. 

Lines  of  constant  maximum  principal  stress  and  strain  in  the  vicinity 
of  the  flap  termination  are  shown  in  Figures  43  and  44,  respectively. 

The  stress  values  a-e  normalized  in  terms  of  <Tr7)  and  the  strain  values 

E  (  rZ  avg 

ri)c\x 

are  shown  in  terms  of  ~r: — : — '■  (a  bond  length  of  94.0  in,  and  a  radial 

x  *  v7  / 

r*-  avg 

distance  of  15.895  in.  were  used  in  calculatin'-  the  average  shear  stress). 
The  stress  values  indicate  a  high  stress  concentration  at  the  flap-insulator 
bond  termination  as  expected.  The  overall  gradient  is  similar  to  that  for 
thermal  shrinkage  loading.  This  similarity  holds  quantitatively,  as  shown 
in  Figure  45,  wherein  the  ratio  of  the  maximum  principal  stress  inside  of 
the  liner/propellant  to  that  at  the  liner-flap  bondline  (H^)  decreases 
almost  identical  to  that  for  shrinkage  loading.  Lines  of  constant  maximum 
principal  strain  are  essentially  the  same  for  axial  acceleration  and  thermal 
shrinkage  loading,  as  shown  in  Figures  7  and  44. 
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Finite-element  stress  solutions  were  obtained  for  W/b  =  0.5,  0.67, 
and  0.8  for  the  cylinder  with'L/D  =  3.D.  Computed  results  for  the  critical 
stresses  at  the  flap  termination  are  shown  in  Figure  ^ 6 .  The  peak  stresses, 
when  normalized  with  respect  to  the  average  case-bond =shear  stress,  are 
essentially  independent  of  web  'fraction.  The  peak  value  of  the  maximum 
principal  stress  (for  a  unit  vulue  of  H<y)  is  approximately  equal  to 
4.6  times  the  average  bond  shear  stress;  whereas,  peak  value  for  the 
radial  stress  is  approximately  4,4  times  the  average  bond  shear  stress. 

The  peak  value  of  the  bond  shear  stress  is  approximately  1.4  times  the 
average  bond  shear  stress.  The  shear  stress  at  the  location  of  the  peak 
radial  stress  is  only  0.17  times  as  large  as  the  peak  radial  stress. 

i  . 

\ 

Finite-element  stress  solutions  for  W/b  =  0.8  were  obtained  for 
L/D  =  1.0,  3.0,  and  5.0.  The  key  stresses  at  the  flap  termination  are 
shown  in  Figure  47  as  a  function  of  1,/D.  These  critical  (normalized) 
stress  parameters  are  essentially  independent  of  length- to-diameter  ratio, 
as  well  as  web  fraction.  Values  for  L/D  =  1.0  are  only  7  percent  higher 
than  those  for  L/D  =3.0  and  5.0.  One  would  anticipate  higher  stresses 
relative  to  the  average  for  shorter  cylinders;  however,'  dependence  upon 
L/D  is  very  moderate  for  values  1.0  and  greater. 

It  appears  that  the  normalized  bond  termination  stresses  for  axial 
acceleration  loading  a:re  basically  dependent  upon  the  local  flap  termina¬ 
tion  geometry  and  are  indeperident  of  L/D  and  W/b!,  In  this  regard,  the 
cylinder  behavior  is  much  different  than  for  thermal  shrinkage  and 
internal  pressure  loadings.  The  peak  bona  stresses  for  axial  acceleration 
are  sensitive  to  i lap-liner  thickness  in  the  same  quantitative  manner  as 
for  thermal  shrinkage  and  internal  pressure  loadings.  No  study  was  per¬ 
formed  for  variations  in  flap  length  or  liner-flap-insulator  modulus. 
However,  since  the  same  stress-strain  patterns  exist  near  the  flap 
termination  for  all  three  loading  conditions,  it  is  expected  that  variations 
for  thermal  shrinkage  loading  are  representative  of  axial  acceleration 
loading  as  well. 
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SECTION  VIII 


SUMMARY  AND  CONCLUSIONS 

The  important  results  and  conclusions  obtained  from  the  parameter 
study' are  as  follows:  1 

(1)  Case-bond  stresses  for  thermal  shrinkage  and  internal  pres¬ 
sure  loading  become  larger  as  web  fraction  and  length-to- 
diameter  ratio  increase.  However,  as  L/D  and  W/b  increase, 
the  end  termination  stresses  decrease  relative  to  the 
centerport  hoop  stress  and  the  radial  bond  stress  at  thd 
mid-cylinder  location. 

(2)  Case-bond  stress  distributions  at  flop  terminations  are  a 

function  of  the  local  geometry  only  and  are  relatively 
insensitive  to  changes  in  W/b  and  L/D  for  the  loadings 
studied.  , 

(3)  Modest  changes  in  flap  and  case-bond  liner  stiffness  have 
little  effect  on  local  case-bond  stress  levels. 

(4)  Flap  length  is  a  less  significant  factor  than  flap  thick¬ 
ness  in  determining  stresses  at  the  flap  terminations. 

(5)  Case-bor.e  and  propellant  stresses  can  be  predicted,  quite 
accurately  for  internal  pressure  loading  using  thermal 
shrinkage  solutions ; 

,  (6)  The  case-b'ond  radial  stress  at  flap  terminations  is  much 
larger  than  the  shear  stress  under  axial  acceleration 
loading,  as  well  as  under  thermal  shrinkage  and  internal 
pressure  loading. 

i  1  ■ 

(7)  Flap  termination  stresses  relative  to  the  average  case- 

bond  shear  stress  are  insensitive  to  L/D  and  W/b  for  axial 
acceleration  loading. 
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Figure  1.  Cylinder  Model  for  Flap  Termination  Studies 


Figure  2.  Finite  Element  Model  for  Cylinder  with  L/D  =1.0  and  W/b 


igure  3.  Finite  Element  Grid  Detail  Adjacent  to  Flap  Termination 


«I!!ATION 


Figure  4.  Case  Bond  Liner  Stresses  for  a  Cylinder  with  L/D  =1.0  and 
W/b  =  0.8  Under  Thermal  Shrinkage  Loading 


FLAP  TERMINATION 


Figure  5.  Case  Bond  Liner  Stresses  Near  the  Flap  Termination  for  a 
Cylinder  with  L/D  =  1.0  and  W/b  =  0.8  Under  Thermal 
Shrinkage  Loading 
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Thermal  Shrinkage  Loading 


Figure  8.  Case  Bond  Liner  Stresses  Near  the  Flap  Termination  for  a 

Cylinder  with  L/D  =1.0  and  W/b  =  0.8  (Ef^gp  =  200  psi)  Under 
Thermal  Shrinkage  Loading 


Figure  11.  Lines  of  Constant  Maximum  Principal  Strain  ^ — g— J  in  the  vicinity 

of  the  Flap  Termination  for  a  Cylinder  with  L/D  =  1.0  and  W/b  =  0 
(Efiap  =  200)  Under  Thermal  Shrinkage  Loading 
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Figure  13.  Propellant  Maximum  Principal  Stresses  at  the  Propel lant- to 
Adjacent  to  the  Flap  Termination  for  a  Cylinder  with  L/D  = 
W/b  =  0.8  Under  Thermal  Shrinkage  Loading 


Liner  Bond 
1 . 0  and 


A3 


Figure  15.  Propellant  Maximum  Principal  Strains  at  the  Propel lant- to-Liner  Bond 
Adjacent  to  the  Flap  Termination  for  a  Cylinder  with  L/D  =1.0  and 
W/b  =  0.8  Under  Thermal  Shrinkage  Loading 


Figure  20.  Ratio  of  the  Maximum  Principal  Stress  in  the  Flap,  Liner,  and 
Propellant  to  that  at  the  Flap-Liner  Interface  in  the  Vicinity 
of  the  Flap  Termination  for  a  Cylinder  with  L/D  =1.0  and 
W/b  =  0.8  Under  Thermal  Shrinkage  Lo'ding 
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Figure  22.  Shear  Stress  in  the  Case  Bond  Liner  for  a  Cylinder  with  L/D  =  1.0 

and  W/b  =  0.8  and  a  Variable  Flap  Length  Under  Thermal  Shrinkage 
Loading 
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Figure  24.  Strain  Energy  in  Cylinder  under  Themal  Shrinkage  Loading,  with 
L/D  =  1.0  and  W/b  =  0.8,  as  a  Function  of  Flap  Length 


Figure  25.  Gradient  in  Maximum  Principal  Stress  as  a  Function  of  Normalized 
Radial  Distance  from  Flap-Insulator  Bond  Discontinuity  for 
Cylinder  with  I./D  =1.0  and  W/b  =  0.8 
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Figure  33.  Case  Bond  Liner  Stress  in  Radial  Direction  for  Cylinder  with 
L/D  =  1.0  and  3.0,  W/b  =  0.8,  Under  Thermal  Shrinkage  Loading 
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ure  35.  Ratio  of  Radial  Liner  Stress  in  Saddle  to  Mid-Cylinder  Va] 
Cylinder  Under  Thermal  Shrinkage  Loading 


Figure  36.  Ratio  of  Liner  Hoop  to  Radial  Stress  at  Mid-Axial  Position  in 
Cylinder  Under  Thermal  Shrinkage  Loading 
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Figure  39.  Case  Bond  Liner  Stresses  Near  the  Flap  Termination  for  a  Cylinder 

with  L/D  »  1.0  and  W/b  ■  0.8  Under  1.0  psi  Internal  Pressure  Loading 


Figure  42.  Case  Bond  Liner  Stresses  Adjacent  to  Flap  Termination  in  a  32-Inch  Inameter 
Cylinder  with  L/D  -  3.0  and  W/b  -  0.8  Under  1  g  Axial  Acceleration 
Loading  (/*w  =  0.064  lb/ in. 3) 
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;ura  43.  Lines  of  Cmstant  Maximum  Principal  Strain  {  — — -  jin  the  Vicinity  of  the  Flap 

\Tavg/ 

Termination  for  a  Cylinder  with  L/D  =3.0  and  W/b  =  0.S  llndar  Axial 
Acceleration  Leading 


AVG 


Figure  ->7 .  Normalized  Case  Bond  Liner  Stresses  Adjacent  to  Flap  Termination  for  Cylinder 
with  T.v/b  =  0.8  Subjected  to  Axial  Acceleration  Loading 


